The three-dimensional Anderson model of localization with binary random potential 
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We study the three-dimensional two-band Anderson model of localization and compare our results 
to experimental results for amorphous metallic alloys (AMA). Using the transfer-matrix method, 
we identify and characterize the metal-insulator transitions as functions of Fermi level position, 
band broadening due to disorder and concentration of alloy composition. The appropriate phase 
diagrams of regions of extended and localized electronic states are studied and qualitative agreement 
with AMA such as Ti-Ni and Ti-Cu metallic glasses is found. We estimate the critical exponents 
I'w, ve and when either disorder W , energy E or concentration x is varied, respectively. All our 
results are compatible with the universal value v ~ 1.6 obtained in the single-band Anderson model. 

PACS numbers: 71.23.An, 71.30.+h, 71.55.Jv 



I. INTRODUCTION 

Amorphous metallic alloys (AMA) offer the possibil- 
ity of continuously changing their composition while at 
the same time avoiding structural phase transformations. 
Thus they allow for systematic studies of their physi- 
cal properties within a single phase as temperature and 
other external control parameters, e.g., pressure, are var- 
ied. Many such investigations (see, e.g., jj) have been 
devoted to the investigation of the electrophysical prop- 
erties (in particular, electrical resistivity) of AMA; these 
studies have revealed considerable differences in the be- 
havior of the electrical resistivity g of AMA in compari- 
son with that of their crystalline counterparts. 

The transport properties of transition-metal-based 
AMA are of special interest because of unique physical 
anomalies, unexplained by conventional transport theory, 
that exist in these materials. Among them are the nega- 
tive temperature coefHcient of resistivity, the Mooij cor- 
relation of resistivity and thermoelectric power — the 
higher the resistivity the lower is the temperature coeffi- 
cient of the resistivity — , the resistivity saturation, the 
sign reversal of the Hall coefficient, a negative magneto- 
resistance, and the breakdown of Mattiessen's rule. For a 
review, we refer the reader to ref. The first attempt 
for a theoretical understanding of the electronic trans- 
port in noncrystalline materials was the Ziman theory 
and its extensions Ij. However, the applicability of the 
Ziman diffraction model for strongly disordered and thus 
high-resistivity amorphous alloys is questionable 0, 1^. 
In addition, the Ziman theory cannot explain, e.g., that 
the changes of the resistivity (giuquid - 6'soiid)/£'soiid at 
the melting point are only 0.01 and 0.09 for Fc and Co, 
respectively Thus it appears that the resistivity of 
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these two elements depends only weakly on their atomic 
structure factor S{k). 

At low temperature T, an even more significant differ- 
ence between the behavior of crystals on the one hand 
and disordered solids on the other hand is seen: suffi- 
ciently strong disorder can give rise to a transition of the 
transport properties from conducting behavior with re- 
sistance i? > to insulating behavior with i? = cxd as 
T — > as was pointed out by Anderson in 1958 This 
phenomena is called the disorder-driven metal-insulator 
transition (MIT) and it is characteristic to non- 

crystalline solids. The mechanism underlying this MIT 
was attributed by Anderson not to a finite gap in the en- 
ergy spectrum which is responsible for an MIT in band 
gap or Mott insulators 0. Rather, he argued that the 
disorder will lead to interference of the electronic wave 
function ip{r) with itself such that it is no longer extended 
over the whole solid but is instead confined to a small part 
of the solid. This localization effect excludes the possibil- 
ity of diffusion at T = so that the system is an insulator. 
A highly successful theoretical approach to this disorder- 
induced MIT was put forward in 1979 by Abrahams et 
al. [lOj . This "scaling hypothesis of localization" details 
the existence of an MIT for non-interacting electrons in 
three-dimensional disordered systems at zero magnetic 
field and in the absence of spin-orbit coupling. 

These ideas have also been applied to the analysis of 
AMA transport properties j^. They provide a back- 
ground for an explanation of the Mooij correlation 
and were successful in the description of the transport 
properties of amorphous semiconductors [T^ . Many pa- 
pers have since discussed the infiuence of localization (or 
quantum interference effects) on the transport properties 
of AMA 11 in in 113. It was shown ^ that Anderson 
localization is responsible for regions of high resistivity. 
The importance of quantum corrections in analyzing the 
conductivity of so-called highly resistive {g > ISO/ificm) 
metallic glasses was demonstrated There are two 

sources of these corrections: the disorder-induced "lo- 
calization effect" and the electronic "interaction effect" . 
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The former caused the observed higher temperature re- 
sistivity. Furthermore, it has been argued [23, |^ that 
for amorphous metals locahzation effects are vahd even 
at room temperature. 

In a recent study [23 the influence of disorder on the 
transport properties of Ti-Ni and Ti-Cu metallic glasses 
was analyzed within the Anderson model of localization. 
The results elucidated qualitatively the correlation be- 
tween the observed details of the electronic structure and 
the behavior of resistivity vs. temperature for this binary 
AMA. It was shown that the temperature cocfhcient of 
resistivity depends on the position of the Fermi level E-p 
and on the (assumed) position of the mobility edge Ec- 
The latter separates extended (conducting) states from 
localized (insulating) states 0, l2,l| | . Unfortunately, the 
value of Ec had to be inferred from the hypothesis that 
the two subbands of the valence band of the binary alloy 
have tails with localized electronic states like the tails of 
the usual one-band Anderson model 6] . Evidently, a cal- 
culation of the position of the mobility edge E^ should be 
carried out for the case of two subbands. This is what we 
intend to do in the present manuscript. In particular, we 
will investigate whether localized electronic states exist 
in the central region between the two subbands as as- 
sumed in [23|. Furthermore, we show how the concentra- 
tion parameter, absent in the usual single-band Anderson 
model, influences the position of the mobility edges. 



II. NUMERICAL APPROACH 

A. The binary Anderson model of localization 

According to the photo-emission data and theoreti- 
cal estimates the valence band of a binary, 
transition-metal-based amorphous alloy can be assumed 
to be a superposition of two valence bands A and B. In 
order to model this within the Anderson model, we use 
the standard Anderson Hamiltonian 

H = J2e,\^){^\+J2u,\^){J\ (1) 

with orthonormal states \i) corresponding to electrons 
located at sites i = {x, y, z) of a regular cubic lattice 
with periodic boundary conditions. The hopping inte- 
grals tij are non-zero only for nearest neighbors and we 
set the energy scale by choosing t.y — t = 1. The two- 
subband structure observed in the experiments is incor- 
porated into the potential energies Si. The randomness 
is modeled by 

1. independent random variations £^(^1) G [ea ~ 
WA/2,eA + Wa/2] and e,(B) e [eb - Wb/^^Eb + 
Wb/2], 

2. random spatial distribution of potential energies ea 
and Eb within the cubic lattice. 



The mean values ea and Eb are chosen according to the 
central energies of the two subbands. The parameters 
Wa and Wb specify the disorder strength in each energy 
band and the random variation of A and B sites models 
the compositional disorder of the binary alloy. Note that 
with this choice of parameters, the band edges of the 
subband A in the limit of large system size are given by 
ea — Q — Wa/2, -I- 6 -I- Wa/2 and similarly for the other 
subband with A ^ B. 

In summary, the model is described by the following 
parameters 

1. eab = Eb ^ Ea, the distance between the centers 
of the two bands (we choose > so that the A 
band is the lower one), 

2. XA, the concentration of A sites (of course, then 

= 1 — XA is the concentration of the B sites) , 

3. Wa and Wb, the widths of the box distributions of 
the potential energy. 

B. Transfer-matrix method 

Since we are interested in the position of the mobil- 
ity edges and thus the localization lengths, we use the 
transfer-matrix method (TMM) jS, iM] for the in- 
vestigation of model The localization length A de- 
scribes the exponential decay of the wave function and 
we compute it using TMM for quasi- ID bars of cross sec- 
tion M X M and length L ^ M. As is customary in 
the TMM approach, the stationary Schrodinger equation 
H^* = E^ is rewritten in the recursive form 

(2) 

^i. Hi, and are wave function, Hamiltonian matrix, 
and transfer matrix of the ith slice of the bar, respec- 
tively. Unit and zero matrices are denoted by 1 and 0. 
Given an initial condition (^j) equation Q allows a re- 
cursive construction of the wave function in the bar geom- 
etry by adding more and more slices. A(Af, w) is then ob- 
tained from the smallest Lyapunov exponent of the prod- 
uct TlTl-1 • • • T2T1 of transfer matrices 0], where the 
length L of the bar is increased until the desired accu- 
racy of A is achieved. With increasing cross section of the 
bar the reduced localization length A]\j{w) = X{M, w)/M 
decreases for localized states and increases for extended 
states. Thus it is possible to determine the critical pa- 
rameter Wc at which Am is constant as a function of the 
varied parameter w — e.g., E, eab, xa, Wa, or Wb — 
from plots of Am versus M . 

C. Finite-size scaling 

The MIT in the Anderson model of localization is ex- 
pected to be a second-order phase transition JQ] . It is 
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characterized by a divergent correlation length ^oo{w) = 
C\w — Wn\~'^, where v is the critical exponent, C is a 
constant Q , w is any of the external control parameters 
given above and Wc is its critical value at which the MIT 
occurs. To construct the correlation length of the infinite 
system ^oo from finite-size data Am [li HE HE H^i the 
one-parameter scaling hypothesis is employed, 



A 



M 



(3) 



All values of Am{w) are expected to collapse onto a single 
scaling curve /, when the system size is scaled by ^oo(w). 
In a system with MIT such a scaling curve consists of 
two branches corresponding to the localized and the ex- 
tended phase. One might determine v from fitting by 
a finite-size scaling (FSS) procedure |23|. But a higher 
accuracy can be achieved by fitting directly the raw data 
|27| . We use fit functions which include two kinds of 
corrections to scaling: (i) nonlinearities of the w depen- 
dence of the scaling variable and (ii) an irrelevant scaling 
variable with exponent —y. Specifically, we fit 



Am = loiXrM'/n + M-yf\{xrM'^n 



(4) 



fn = y^Q»i 



Xr(t^) = W 



E 

n=2 



(5) 



with Lo — \wc ~ w\/wc and expansion coefficients a„i and 
bn- Choosing the orders and m-^ of the expansions 
larger than one, terms with higher order than linear in 
the w dependence appear. This allows to fit a wider w 
range around Wc than with the previously used linear 
fitting |3lj| . The linear region is usually very small. The 
second term in equation ^ describes the systematic shift 
of the crossing point of the Am{w) curves [30ll3ll|. 

In the present case of the two-band Anderson model, 
we can in principle have an MIT as a function of 

1. energy E for fixed Wa, Wb, £ab and xa, 

2. disorder strengths Wa, Wb for fixed E, eab and 

3. concentration xa for fixed E, sab and Wa, Wb- 

Due to universality, we expect the corresponding critical 
exponents ve, ^Wai ^Wb and to be the same. Addi- 
tionally for each control parameter we test whether the 
fitted values of Wc and v are compatible when using dif- 
ferent expansions of the fit function, i.e., different orders 
rir and rrir (3Z |. 



III. RESULTS AND DISCUSSION 

For the determination of the mobility edge (2^, we 
have determined the localization lengths from TMM with 
an error equal to or less than 1% for cross-sections up to 
M = 12. For the computation of the critical exponents. 
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FIG. 1: Dependence of the reduced localization length Am 
on energy E for different system sizes M = 6, 8, 10 and 12 at 
EAB = 17, XA = 0.5 and Wa = Wb ~ 5. The two shaded 
regions indicate extended states, the dashed lines the position 
of the 4 mobility edges, and the vertical solid lines the band 
centres of the subbands. No error bars are shown because 
the accuracy of the Ajvj values is better than the size of the 
symbols. 



we have used the available 1% data and generated ad- 
ditional data with 0.05% error up to at most 0.1% at 
selected points in the phase diagrams for high-precision 
estimates. 

In Fig.^ we show a typical dependence of the reduced 
localization length Am on energy for different system 
sizes M. The two peaks correspond roughly to the po- 
sitions of the two subbands at ea and Eb, respectively. 
This is in accordance with the theoretical and experi- 
mental results [23, |3| mentioned in the introduction as 
motivation of our studies. Note that for the chosen pa- 
rameters, the density-of-states limit given in section Hi Al 
for the highest energy state in the lower subband coin- 
cides with the smallest possible energy state in the upper 
subband at ea + Q + Wa/^ = es - 6 - Wb/2. At the 
positions indicated by the dashed lines in Fig.^ we see a 
reversal in the systematic size dependence of the reduced 
localization lengths. Therefore these indicate the mobil- 
ity edges and we observe localized electronic states at the 
band edges and in the central region between the A and 
B subbands. We remark that the Am{E) curves demon- 
strate a peculiar feature: the positions of the maximum 
values of Am do not coincide with the central energies 
ea and es of the two subbands. Thus although the den- 
sity of states between the two bands is appreciable, we 
observe the surprising fact that the region of extended 
states in each band has shrunk. 

In Fig. 121 we show high precision data at the upper 
mobility edge of the A subband at £'c « 2. Again the 
accuracy of the data is so high that error bars would 
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FIG. 2; FSS plot of the reduced localization length Am near 
the mobility edge i5c ~ 2 of Fig. The solid lines are fits of 
the data according to equations I^J and @ with rir = 1 and 
rur = 1. The inset shows the scaling function corresponding to 
the fit and the Am data with symbols o, ■, A, <, V, >, +, *, •, □ 
denoting energies 1.9, 1.92, 1.94, . . . , 2.08, 2.10, respectively. 



localized electronic states in the central region can be 
found. In Fig. ^ we show the mobility edges for differ- 
ent eab- We note that the mobility edge for the usual 
Anderson model equals approximately ±6.2 at disorder 
= 4 H, H^l ■ This agrees with the result for eab = 
in Fig. H 



C. Energy-concentration phase diagram 

Obviously, the values for Am and the mobility edges 
should strongly depend on the concentration. This is 
indeed the case as we show in Fig. |5| Already a 10% 
difference in atomic composition leads to a pronounced 
asymmetry of the Am curves (compare, e.g., with Fig.^l. 
The energy-concentration phase diagram for eab = 8, 
and Wa = Wb = 4 is presented in Fig. El For xa < 
0.32 {xb > 0.68), all electronic states of subband A [B) 
are localized. Figs. |S1 and |H| correspond to the case of 
overlapping bands such that the center of subband A 
coincides with the edge of subband B and vice versa. 



be smaller than the symbols. Data and FSS analysis 
of similar quality will be used when estimating critical 
exponents in the following. 

A. Energy-disorder phase diagram 

Let us now investigate how the positions of the mo- 
bility edges change when Wa and Wb are changed. We 
set Eab — 17 and xa — xb — 0.5. For convenience we 
choose Wa = Wb- This leads to a symmetry for the 
energy dependences of Am{E) between lower and upper 
subband as shown in Fig. Consequently, the phase 
diagram is symmetric with respect to ea + £ab/'^- We 
find extended electronic states in the vicinity of ea and 
for Wa <Wc'^ 6.5 (instead of 16.5 as for the usual An- 
derson model) as presented in Fig. 13 For larger Wa or 
larger \E — ea\, the states are localized. 

Next, we keep Wa = 4 fixed and vary Wb from to 
8. We find for eab — 17 that the position of mobility 
edges in the subband A is only slightly modified. The 
same behavior is observed for eab — 8, when the overlap 
of the two subbands is large. 

B. Energy-£As phase diagram 

Keeping xa = 0.5 and fixed disorder Wa = Wb = 4, 
we now vary the subband spacing eab- From an ex- 
perimental point of view, this is the difference between 
the atomic constituents of the alloy or, more precisely, 
between their ionization energies. For a variety of eab 
values we have obtained Am curves which are qualita- 
tively very similar to Fig. ^ However, for eab < 6, no 



D. Concentration-eAs phase diagram 

The concentration-e^iB phase diagram is presented in 
Fig. □ for = = and Wa ^ Wb = 4. We see 
that the critical concentration is larger than the classical 
percolation threshold xa = 0.3116 for large eab- We in- 
terpret this as being due to the localization effect in this 
quantum situation: a larger cluster of connected sites is 
needed to support extended states on sites of the sublat- 
tice A (or B). For very large eab = 100 and 10'^ (not 
shown) the critical concentration decreases to 0.416 and 
0.384, respectively, but remains larger than in the classi- 
cal case. For large subband spacing sites with onsite po- 
tential E B can be viewed as potential barriers to electrons 
from subband A and vice versa. Thus electrons with en- 
ergy « Ea have a low probability to hop onto B sites 
and we are effectively studying the case of percolation. 
But even though eab might be large, there will always 
remain a finite but small probability for tunneling thus 
distinguishing our model from classical percolation. And 
since we are studying the wave equation , localization 
effects not present in classical percolation are important 
and lead to insulating behavior even when the concen- 
tration of sites supports a classically percolating cluster. 
Therefore, we find that the threshold for transport (ex- 
tended states) in the quantum case is higher than for the 
classical one. We note that similar effects have been ob- 
served previously in studies of quantum percolation |34j . 

On the other hand, for small EaBi the critical concen- 
tration becomes less than the percolation threshold since 
tunneling between A and B sites is more significant. The 
phase diagram in Fig. |7| does not exactly show the mo- 
bility edges at which all extended states disappear for 
a given eaBi because it is determined at E — 0. But 
as shown in Figs. ^ and El the A band is not symmet- 
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FIG. 3: The energy-disorder phase diagram for eab = 17 and xa = 0.5. The shaded region indicates extended states and the 
horizontal line denotes the value of ea- Here and in the following phase diagrams, □ and o indicate lower and upper mobility 
edges for the A band, and, if included, upper and lower mobility edges for the B band. 
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FIG. 4: The energy-SAs phase diagram for Wa = Wb = 4 
and Xa = 0.5. The shaded region indicates extended states. 
The upper (B) and lower (A) bands are symmetric with re- 
spect to eab/^ and the B band has been constructed using 
this symmetry. 
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FIG. 5: The energy dependence of Am for a majority of A 
sites (concentration xa ~ 0.55). System sizes are M — 6,8, 10 
and 12, eab ~ 8, and Wa = Wb ~ 4. The two shaded regions 
indicate extended states. The solid and dashed vertical lines 
denote the band energies ea and eb and the position of the 4 
mobility edges, Ec, respectively. 



ric with respect to eyi(= 0). As a consequence, more 
extended states occur for energies just below ea- This 
means that the phase boundary in Fig. [T] will be shifted 
slightly downwards, if E is varied, too. 

We note that an analogous phase diagram has been 
obtained recently [s^ for the case of a two-dimensional 
binary alloy. The shape of the curve separating local- 
ized and extended electronic states is quite similar to our 
result shown in Fig. \7\ 



E. Disorder-concentration phase diagram 

The concentration dependence of the critical disorder 
— the disorder at which all states for a given energy 
become localized — is presented in Fig. |S1 for energy 
E = Ea- One can see that the critical disorder strongly 
depends on the concentration. For xa — 1, we recover 
the result Wc = 16.5 of the single-band Anderson 
model. No extended electronic states are observed for 
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FIG. 6: The energy-concentration phase diagram at band 
separation eab = 8 and disorders Wa = Wb ~ 4. The 
shaded regions denote extended states. The upper region of 
extended states has been constructed from point symmetry 
with respect to eab /2 and xa ~ 0.5. The thick solid lines and 
the triangle indicate the inferred concentration dependencies 
of 4 different AM A as discussed further in section liVI 
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FIG. 7: The concentration-£As phase diagram ai E = ea ~ 
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percolation threshold value xa = 0.3116. 



15 



10 



CD 
T3 

O 
CO 



localized 



extended 



0.2 0.4 0.6 

concentration (x^) 



0.8 



FIG. 8: The disorder-concentration phase diagram for eab = 
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extended states. The arrow refers to the discussion in the text 
as well as to the data presented in Fig. |^ The x denotes the 
single-band result. 
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FIG. 9: The disorder dependence of Am for system sizes M = 
6, 8 and 10 at eab = 17, Wis = 4, £ = and xa = 0.38. The 
shaded region indicates extended states and the two vertical 
lines denote approximate estimates of critical disorders Wc- 



Xa < 0.35 in agreement with the idea that for small con- 
centration the electronic states of these atoms form the 
usual localized donor subband of separated impurities. 

For the concentration interval from 0.35 to 0.40 we 
obtain an additional transition as a function of disorder 
at small Wa besides the usual Anderson MIT. For exam- 
ple, as shown in Fig.|51for xa = 0.38 we observe localized 
behavior of Am{M) for small disorders. Increasing the 
disorder beyond Wa > 0.7, we see extended states. At 
Wa ^ 3.2, the character of the states changes back to 



localized. This behavior can be understood as follows: 
for Xa ~ 0.38 as indicated by the arrow in Fig. |S1 the 
system at Wa = represents a binary alloy with large 
energy separation eab = 17. All states are localized, 
i.e., non-quantum percolating, on an already classically 
percolating set of sites {i(^)} in the A band. Small ad- 
ditional potential disorder will lead to an increase of cer- 
tain ei(^A) as well as a decrease in certain £hb)- A small 
change of the scattering conditions in these sites can then 
quickly allow enhanced transport across a now quantum- 
percolating backbone. But upon further increasing the 
disorder, the localization in each band will quickly lead 
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will of course be purely qualitative. Nevertheless, it al- 
ready suffices to understand much of the transport prop- 
erties of such AMA. E.g., as mentioned in section it 
is the position of the mobility edge relative to the value 
of the Fermi energy that determines the properties of 
transition-metal-based AMA. In this spirit, the Mooji 
correlation is due to Ep being situated in a region of 
localized states, whereas the weak changes in resistivity 
at melting for Fe and Co as well as the positive temper- 
ature coefficient of resistivity for Fe-based AMA result 
from Ep being located in the center of the 3d-electronic 
band and hence, in the region of extended states. 
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FIG. 10: Concentration dependence of the critical exponent 
vw obtained when varying the disorder across the mobility 
edge of Fig.lSl i.e., for E = ea = 0, £ab = 17 and Wb = 4. 
The different symbols represent various orders rir and rrir (in- 
dicated in the legend) of the expansion (|1J used for the fitting. 
The error bars denote one standard deviation as obtained 
from the non-linear fitting procedure. 



to localization again. 



F. Value of the critical exponent 

The critical exponent for the present model can be esti- 
mated by crossing the mobility edges as function of either 
disorder, energy or concentration, giving rise to the three 
exponents and v^^ respectively. Due to universal- 

ity, their values should coincide with the critical exponent 
J/ « 1.6 of the single-band Anderson model since the uni- 
versality class remains unchanged by the introduction of 
the binary disorder 36] . 

We have estimated the critical exponent v for almost 
all points in the phase diagrams 13 E El El and |S1 as de- 
scribed in section III CI For example, we compute v\\r 
along the mobility edge of the disorder-concentration 
phase diagram of Fig. |S| In Fig. ^| we show the re- 
sulting estimates when using different orders and 
for the fit function Q). As usual, the spread in val- 
ues is somewhat larger than the least-square-error bars 
seem to suggest [s^- Summarizing the results, we find 
vw = 1.56 ± 0.06, VE = 1-64 ± 0.05 and = 1.60 ± 0.07, 
compatible with « 1.6 [3^. 



IV. DISCUSSION AND COMPARISON WITH 
EXPERIMENTS 

Let us now compare the above results of the binary An- 
derson model with experimental measurements of trans- 
port properties of AMA [2^ Is^ . Such a comparison 



A. Transport properties of Nii-^jPi, 

The ionization energy of P is higher than that of Ni, 
so our A band corresponds to P atoms and the B band 
to Ni atoms with eb > £a as, e.g., shown in Fig. 
According to the theoretical and experimental results of 
Ref. Ey of the amorphous alloy Nii-xPa; is located 
near the upper edge of the Ni-3d band (the B band in our 
study). When the P concentration xp increases from 0.15 
to 0.27 the temperature coefficient of resistivity (TCR) 
decreases from positive to negative values ,38] . The MIT 
is observed near xp = 0.25. In Fig. we indicated the 
corresponding region in the phase diagram. The two- 
band, binary random potential Anderson model can in- 
deed qualitatively capture the observed experimental re- 
sults. Even a quantitative comparison might be possible 
if a reliable way of obtaining values for ea, £b, and Wa 
and Wb becomes known. Note, however, that in the 
experiments, it is the hopping between A and B atoms 
at their respective random positions that leads to the 
broadening of the A and B bands and we have modeled 
this only qualitatively in the present work by the onsite 
disorders Wa and Wb- 



B. Transport properties of Tii-aj-yCua, Ay and 

Til — X — y N^lx A.y 

A previous experimental study of the electronic struc- 
ture of Tii_a;_yCua: Aj, and Tii_2;_j,Nia: Aj, amorphous al- 
loys [23| (the symbol Aj, represents additional admix- 
tures) has shown that the valence band has two main 
peaks. The lower peak is formed predominantly by the 
3d states of Cu or Ni whereas the upper peak is due to 
the 3d states of Ti. The Fermi level is located in the cen- 
tral region between the peaks and is hardly shifted upon 
changing the concentrations. 

Within the present two-band model, this experimental 
situation may be modeled as shown in Fig. El The line for 
the Tii_a;_j,Nij;Aj, alloys is situated somewhat below the 
one for Tii^x-yCuxAy because Ni has one electron less 
than Cu. We believe that it is this difference that leads to 
the smaller value of the negative TCR for Tii_2;_yNia;Aj, 
when compared to Tii-x-yCuxAy. In both cases, E-p 



8 



TABLE I: Temperature coefficient of resistivity Q300 at 300 K 
for Tii_j;_yCuj; Ay and Tii-^-yNi^Ay amorphous alloys |2^ . 



AMA 


aaoo 


AMA 






iq-Vk 




10-*/K 


Ti62CU33P5 


-0.70 


TiroNiasSis 


-2.21 


Ti47CU45Ni5Si3 


-1.11 


Ti6oNi36P2Si2 


-2.72 


Ti48Cu45Ni5P2 


-1.05 


Ti5oNi45P5 


-2.37 


Ti46CU45Ni5Si2P2 


-3.03 


Ti45Ni5oP5 


-3.28 



lies in a region of localized states in agreement with the 
negative TCR. In Table H we show the values obtained 
for TCR. 

Let us now analyze the concentration dependence of 
the TCR. When the Ti concentration decreases, the TCR 
also decreases as shown in Table This behavior can be 
explained from the concentration-energy phase diagram 
of Fig. El by the increase of the mobility edge . This 
in turn leads to an increasing distance between and 
the region of unoccupied extended electronic states. The 
states close to Ep become more localized, leading to a 
more negative TCR. 

Therefore, at least at a qualitative level, the present 
two-band Anderson model explains not only the nega- 
tive TCR of Tii„j;_j^Cua;Ay and Tii_2:_j^Ni^Aj^ glasses 
but also the changes in TCR due to changes in com- 
position (Cu or Ni) and concentration. Furthermore, 
metallic glasses with a "metal" -like conductivity can be 
also treated within this model. For example, the Fermi 
level of Fe86Bi4 alloy is situated in the region of extended 
states as indicated in Fig.Elin agreement with its positive 
TCR (23]. 



C. Overlapping the two bands 

Thus far, motivated by the experimental results, we 
concentrated on the case of separated bands, i.e., large 
Sab- To make the two regions of extended states in Fig.|Bl 
overlap, we can in principle (i) keep eab fixed and change 
the disorder, or (ii) keep fixed disorder and decrease £ab 
as in Fig. 21 From Fig. O we know that the position 
of mobility edges does not change a lot when decreas- 
ing Wa- On the other hand, increasing Wa will lead to 
a narrower region of extended states in Fig. and no 
overlap. Therefore, we choose to reduce sab in order to 
study what happens when the two bands overlap in the 
energy-concentration diagram. 

In Fig. 111! we show such an overlap for eab = 6 and 
Wa = Wb = 4. This situation corresponds to large 
overlap of the two subbands. We see that the shape of 
the regions of extended states remains nearly unchanged. 
This leads to non-monotonicity in the behavior of the 
upper and lower mobility edges. We remark that in order 
to resolve this behavior, the accuracy of the TMM had to 




0.2 0.4 0.6 0.8 1 

concentration (x^) 



FIG. 11: The energy-concentration phase diagram at band 
separation eab ~ 6 and disorders Wa = Wb = 4. The 
shaded region denotes extended states. The horizontal lines 
indicate eb, eab/ 2 and eb from top to bottom, respectively. 
The upper mobility edge has been constructed from point 
symmetry with respect to eab/2 and xa = 0.5. 



be increased significantly up to very small errors of 0.05% 
for the localization lengths. We emphasize that a similar 
accuracy was needed to resolve the non-monotonicity in 
the disorder-concentration diagram in Fig. [HI 

When further decreasing eab to 5, the behavior of the 
mobility edges becomes regular again and their values 
eventually reduce to the values of the well-known phase 
diagram of the single-band Anderson model ^3 ■ 



V. CONCLUSIONS 

We have shown that the three-dimensional Anderson 
model of localization with binary random potential disor- 
der allows to explain not only various peculiarities of the 
transition-metal-based amorphous alloys but also how 
these change as the chemical composition changes. Even 
for a comparatively large overlap of the two subbands, a 
region of localized electronic states exists in the central 
part between the two subbands. This confirms the as- 
sumption made in Ref. j22| for the experimental results 
in AMA. Even a quantitative analysis, using suitably ex- 
tracted parameters from density-of-states measurements 
appears possible. Of course in this case ternary and qua- 
ternary AMA will have to be distinguished by additional 
onsite potentials similar to Sa and Sb for the binary 
model. 

In addition to the usefulness of the two-band Ander- 
son model for a qualitative comparison with experiments, 
there is also the interesting fact that the model allows to 
study the concentration dependence and the associated 
MIT. This concentration dependence is clearly a promi- 
nent feature present in most disordered materials. We 
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find that, as expected from universality, the universal- from the DAAD via Ref. 322, A/01/20012 and the DFG 
ity class of this additional transition remains unchanged via SFB393 . 
with respect to a single-band Anderson model. 
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